Critical behavior of interacting two-polymer system 
in a fractal solvent: an exact renormalization group 
approach 

I Zivicf, S Elezovic-HadzicJ and S Milosevic! 

fFaculty of Natural Sciences and Mathematics, University of Kragujevac, 34000 
Kragujevac, Serbia 

^Faculty of Physics, University of Belgrade, P.O.Box 368, 11001 Belgrade, Serbia 
E-mail: ivanz@kg . ac .yu, sukiOf f . bg. ac . yu, emilosevOetf . bg. ac . yu 

Abstract. We study the polymer system consisting of two polymer chains situated 
in a fractal container that belongs to the three-dimensional Sierpinski Gasket (3D SG) 
family of fractals. Each 3D SG fractal has four fractal impenetrable 2D surfaces, which 
are, in fact, 2D SG fractals. The two-polymer system is modelled by two interacting 
self-avoiding walks (SAWs), one of them representing a 3D floating polymer, while 
the other corresponds to a chain adhered to one of the four 2D SG boundaries. We 
assume that the studied system is immersed in a poor solvent inducing the intra- 
chain interactions. For the inter-chain interactions we propose two models: in the 
first model (ASAWs) the SAW chains are mutually avoiding, whereas in the second 
model (CSAWs) chains can cross each other. By applying an exact Renormalization 
Group (RG) method, we establish the relevant phase diagrams for 6 = 2,3 and 6 = 4 
members of the 3D SG fractal family for the model with avoiding SAWs, and for 
b = 2 and 6 = 3 fractals for the model with crossing SAWs. Also, at the appropriate 
transition fixed points we calculate the contact critical exponents, associated with the 
number of contacts between monomers of different chains. Throughout the paper we 
compare results obtained for the two models and discuss the impact of the topology of 
the underlying lattices on emerging phase diagrams. 
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1. Introduction 

The self-avoiding walk (SAW) is well-disposed as a standard lattice model for a flexible 
linear polymer chain in various types of solvents [1]. In this model, the monomers that 
comprise a polymer chain are related to the steps of a random walk that must not 
contain self-intersections, while the surrounding solvent is represented by an underlying 
lattice. In a good solvent, with each step of the SAW we associate the same weight 
factor x, while in a poor solvent, when two non-consecutive monomers of a polymer 
chain become nearest neighbors, we introduce the additional statistical factor u, which 
corresponds to the intra-chain energy e u < 0. Even though an isolated polymer chain 
is difficult to observe experimentally, numerous studies of the single-chain statistics 
have been upheld as a requisite step towards understanding the statistics of many-chain 
systems. A natural extension of a single polymer concept is a model of two interacting 
linear polymers, which may be relevant to perceive behavior of multicomponent polymer 
solutions [2] . To study the critical properties of the two-polymer system we shall apply 
the following two models: The first is the model of two mutually avoiding self-avoiding 
walks (ASAWs), whose paths on a lattice cannot cross each other, and the second is 
the model of two mutually crossing self- avoiding walks (CSAWs), that is, the case of 
the two SAWs whose paths can intersect each other. Various types of models with 
two avoiding SAWs have been successfully applied in the studies of phase transition of 
diblock copolymers [3,4], as well as in the studies of unzipping double-stranded DNA 
molecules [5-8]. On the other hand, the model with two crossing SAWs was applied 
for studying the collapse transition of two-chain interacting system on three- and four- 
simplex lattice [9,10], and Euclidean lattices [11], as well as to study two randomly 
interacting directed polymers on diamond hierarchical lattice [12,13]. 

In this paper we apply both ASAWs and CSAWs model to study the two-polymer 
system that displays both intra- and inter-chain interactions, on the three-dimensional 
(3D) fractal lattices, which belong to the Sierpinski gasket (SG) family of fractals. We 
assume that one of the two polymers is a floating chain in the bulk of a 3D SG fractal, 
while the other is a polymer chain that stays affixed to one of the four boundary surfaces 
(being actually 2D SG fractals) [14]. In the ASAWs model we assume that two SAWs 
are in contact when they approach each other at the distance which is equal to a lattice 
constant, and in this situation we ascribe the contributing contact energy e v to the total 
model energy. Similarly, in the CSAWs model we assume that each crossing between 
two SAW paths corresponds to a contact of two monomers that belong to different 
polymer chains, and therefore we associate the contact energy e c with such a crossing. 
Since in both models, one of two polymers is adhered to one of the fractal boundary 
surfaces, and because its monomers take effect of surface interacting points for the bulk 
floating polymer chain, the proposed models may also be of interest for the problem of 
surface-interacting polymer chain in homogeneous [15-17] and disordered media [18]. 

The main goal of this study is to establish phase diagrams in the space of interaction 
parameters (which consists of the intra- and inter-chain interaction energy parameters), 
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for both models, as well as to calculate the contact critical exponents that describe 
behavior of numbers of monomer- monomer contacts between two polymer chains. 

This paper is organized as follows. In section [2] of the paper, we first describe the 
3D SG fractals for general scaling parameter 6, as well as the ASAWs model. Then, we 
present the general framework of an exact renormalization group method, within the 
model, and elaborate on the phase diagrams, obtained for the fractals designated by 
6 = 2,3 and 6 = 4. We also display our findings for the contact exponents (associated 
with the number of contacts between the two SAWs). The CSAWs model is described 
in section [3j Again, by applying an exact RG method, which is in the latter case 
technically more complicated, phase diagrams and contact exponents for 6 = 2 and 3 
SG fractals are obtained, and discussed. Brief summary and the concomitant conclusion 
are presented in section HI Explicit form of the RG equations for particular fractals are 
given in appendices. 

2. The model of two evading self-avoiding walks 

In this section we are going to apply the renormalization group (RG) method to the 
model of two mutually avoiding self-avoiding walks on the 3D SG family of fractals. 
First, we give a summary of the basic properties of these fractals. We start with recalling 
the fact that each member of the 3D SG fractal family is labeled by an integer 6 > 2 and 
can be constructed in stages. At the first stage (r = 1) of the construction there is a 
tetrahedron of base 6 that contains 6(6+ l)(6 + 2)/6 upward oriented unit tetrahedrons. 
The subsequent fractal stages are constructed recursively, so that the complete self- 
similar fractal lattice can be obtained as the result of an infinite iterative process of 
successive (r — > r + 1) enlarging the fractal structure 6 times, and replacing the smallest 
parts of enlarged structure with the initial (r = 1) structure. Fractal dimension df 
of the 3D SG fractal is equal to df* = ln[6(6 + 1)(6 + 2)/6]/ln6. Each of the four 
boundaries of the 3D SG fractal is itself a 2D SG fractal, with the fractal dimension 
df = ln[6(6+ l)/2]/ln6. 

In the terminology that applies to the SAW, we assign the weight 23 to a step of 
the SAW in the bulk (3D SG fractal), which represents a floating polymer (we mark 
it by P 3 ), and the weight x 2 to a step of the SAW performed on one of the fractal 
boundaries (2D SG fractal), which represents a 2D surface- adhered polymer (marked by 
P2), whose monomers act as interacting counterparts for monomers of the 3D polymer 
chain. To describe the intra-chain interaction of P3 chain, we introduce the Boltzmann 
factor u = e _eil//fesT , where e u < is the interaction energy of two non- consecutive 
neighboring monomers of P3. 

In ASAWs model the two SAWs, that represent polymer chains, must not intersect 
each other. We assume that monomers, belonging to different chains, interact when 
they reach a distance which is equal to a fractal lattice constant, and to a such mutual 
position of P3 and P2 monomers we associate the weight factor v = e~ ev ^ kBT (see figure 
[U^a)), where e v < is the appropriate inter-chain interaction energy. To describe exactly 




Figure 1. The structure of the three-dimensional SG fractal, for b = 2, at the first 
stage of construction, with an example of the bulk polymer chain (P3) depicted by 
green line and the surface-adhered polymer chain (P 2 ) depicted by yellow line. The 
shaded area represents the adhering surface (the two-dimensional SG fractal). The 
intra-chain interactions u, for the P3 polymer, are depicted by blue bonds. In the 
ASAWs model (a) the SAW paths, representing P3 and P2 polymers, cannot intersect 
each other, and two SAWs interact when approach each other at a distance which is 
equal to a lattice constant (red bonds, weighted with v). On the other hand, in the 
case of CSAWs model (b), the polymers P3 and P2 are cross-linked at the two sites, so 
that each contact contributes the weight factor w, while the red bonds (marked by t) 
correspond to the interactions between those monomers which are nearest neighbors 
to the cross-linked points. The two depicted examples for ASAWs (a) and CSAWs 
models (b), contribute the weights x\x2U A v 12 and a^fa^u^ 3 , respectively. 



all possible configurations of the two-chain polymer system within the adopted model, 
we need four restricted partition functions A^'\ B^ r \ and D^ r \ which are defined 

as 

= J2 A {r) (N 3 ,L)x^u L } B {r) = B {r) (N 3 ,L)x^u L } 

N 3 ,L N 3 ,L 

C (r) = J2c {r) (N 2 )x%\ D^= Yl rt r \N 2 ,N 3 ,L,M)x^x^u L v M , (2.1) 

N 2 N 2 ,N 3 ,L,M 

where A^ r \ B^ r \ C^ r \ and T>^ represent the numbers of particular configurations, 
consisting of one or two SAW strands on the r-th fractal structure (see figure [2]). 
For instance, T>( r \N2, N3, L, M) is the number of configurations consisting of A^-step 
P 3 chain with L pairs of non- consecutive nearest-neighbor monomers, and iV 2 -step Pi 
chain, such that there are M contacts between these two chains. The recursive nature of 
the fractal construction implies the following recursion relations for restricted partition 
functions 



A' = 


]T a(N A ,N B )A N *B N », 


(2.2) 




N A ,N B 




B' = 


£ b(N A ,N B )A N *B N *, 


(2.3) 




N A ,N B 




c = 




(2.4) 




N C 




D' = 


£ d(N A , N B , Nc, N d ) A Na B Nb C Nc D n ° , 


(2.5) 




N A ,N B ,N C ,N D 
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Figure 2. Schematic depiction of restricted generating functions used in the 
description of all possible two-SAW configurations, within the r-th stage of the 3D 
SG fractal structure, for ASAWs model. The 3D floating chain is depicted by green 
line, while the 2D surface-adhered chain is depicted by the yellow one. The interior 
details of the r-th stage fractal structure, as well as details of the chains, are not 
shown (for the chains, they are manifested by the wiggles of the SAW paths). The 
functions A^ r \ B^ r \ and C^ r \ describe one-polymer configurations (they are the same 
for both ASAWs and CSAWs models), while the function depicts the inter-chain 
configurations of ASAWs model. 



where we have used the prime symbol as a superscripts for (r + l)-th restricted partition 
functions and no indices for the r-th order partition functions. These relations can be 
considered as the RG equations for the problem under study, with the initial conditions 



which correspond to the unit tetrahedron^. 

Equation (12.41) . alone, describes a single SAW on 2D SG fractal, whereas (12.21) and 
(12. 3D are RG equations for a single SAW on 3D SG fractal. Critical properties of the 
SAW, based on the analysis of these equations, have been well established previously, 
and here we recall their basic properties relevant for the present work. 

First, we describe the behavior of a single 2D SG chain. The RG equation (12.41) . for 
any b, has only one non-trivial fixed point C*, corresponding to the extended polymer 
phase [19,20], that is, the 2D SG chain is always swollen, and it cannot be in the compact 
phase. The corresponding eigenvalue X U2 of (12.41) is larger than 1, and determines the 
value of the critical exponent i> 2 = In bj In X V2 , that governs the behavior of the mean 
end-to-end distance of 2D SG chain (R) ~ (N 2 ) V2 , where (N 2 ) is the average number of 
2D SG SAW steps. 

In what follows we provide short summary of the results concerning the critical 
behavior of a solitary 3D SG chain. Depending on the value of the intra-chain interaction 
parameter u, a single 3D SG chain can be found in three phases: extended chain (for 
u < Ug), ^-chain (when u = Ug) and globule {u > u e ). These phases (for arbitrary b) 
are described by the fixed points (Ae, Be), (Ag, Be) and (Ac, Be), respectively [21-23]. 
The mean end-to-end distance (R) ~ (N 3 ) Ui of SAW on 3D SG fractal, is equal to 
z/ 3 = ln&/lnAj, 3 , where X U3 is the largest eigenvalue of the linearized RG equations 
(12. 21) and (12. 3p . at the corresponding fixed point. For each 3D SG fractal, the following 

X Such initial conditions imply that a SAW can traverse unit tetrahedrons along only one, or two 
nonconsecutive edges. This restriction of the standard SAW model does not alter the critical behavior 
of the system. 



A® = x 3 , B^ 




C<°> = x 2 , D<® = x 3 x 2 v 4 , 



(2-6) 
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relationship > v\ > , is valid, where , u% and u% are the end-to-end distance 
critical exponents in extended, 9 and globule phase, respectively. 

The interacting configurations of P2 and P3 chains are described with the restricted 
partition function D^) . The mean number of contacts between P 2 and P3, on the r-th 
stage of fractal construction, is equal to 

= 4t T, MV^'*> L v m = . (2.7) 

On the other hand, taking into account the function dependance £)( r+1 ) = 
D^ r+1 ){A { r\B { r\C { ~ r \D^), and the fact that A<?\ P (r) , and do not depend on 
the interaction parameter v, we have 

_ 9£)(H-1) g£)(r) 

from which follows that, in the vicinity of the transition fixed point (A*, B* , C* , D*) of 
the two-polymer system, the mean number of contacts (M^), for large r, behaves as 
(MM) ~ \ r D , where 

A -(w) • (2 - 9) 

is relevant eigenvalue of RG equation (I2.5p . calculated at the transition fixed point. 
Knowing that (iV^) ~ X r U3 , one obtains ln(MW)/ln(^ r) ) ~ lnA^/hiA^, i.e. the 
following scaling relation is satisfied 

(MW) ~ (AT 3 (r V , (2.10) 

where 

In Ad 



In A„ 3 



(2.11) 



is so-called contact critical exponent. 

To establish the exact forms of RG equations, for each fractal, one needs to find 
the coefficients a, b, c, and d, that appear in fl2.2l) - fl2.5p . Using the computer facilities, 
by direct enumeration and classification of all possible SAW configurations on the first 
stage of fractal construction, it is feasible to find these coefficients for fractals labelled 



by b = 2, 3 and 4 (see Appendix A). Precise numerical analysis of the obtained RG 



equations (for b = 2, 3, and 4) reveals that two-polymer system can reside in several 
phases, depending on the values of the interaction parameters u and v. In particular, for 
each value of u, there is a critical value v = v c (u), such that for v < v c (u) the two chains 
exist almost independently in the solution. This is indicated by the fact that (A*,B*), 
and C* retain their fixed values that correspond to the solitary chain on 3D SG, and 
2D SG, respectively (see table [1]), and confirmed by calculating the mean number of 
contacts (M^)) between the chains, which quickly approaches some constant value as 
r — > 00. For v = v c (u) fixed values (A*, B*) and C* remain the same as for v < v c (u), 
but D* becomes larger, and (M^) increases with r, obeying the scaling relation (I2.10p . 
with eigenvalue A^ being larger than one. Although there are large number of contacts 
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Table 1. Coordinates of transition fixed points (A*, £>*, C*, D*), obtained via 
renormalization group approach, for ASAWs model on 3D SG fractals labelled by 
b = 2, 3 and 4. Also, we give the corresponding relevant eigenvalues A„ 3 and Xd, 
together with the pertaining values of contact critical exponents <j>. For all values of 
b studied, when v < v c (u), eigenvalue A^ is not relevant, and the mean number of 
contacts between the chains is finite (thus the values of A^ and <ft are n °t given for 
these cases). 



b 


V 


A* 


B* 


C* 


D* 


A„ 3 


Ad 









extended 3D chain (u 


< ue) 










v < v c (u) 


0.4294 


0.0499 


0.6180 


0.1165 


2.7965 


< 1 




2 


v = v c (u) 


0.4294 


0.0499 


0.6180 


2.3303 


2.7965 


2.0904 


0.7170 




v > v c {u) 


0.4294 


0.0499 





3.0887 


2.7965 


2.9537 


1.0532 




v < v c (u) 


0.3420 


0.0239 


0.5511 


0.0779 


5.3620 


< 1 




3 


v = v c (u) 


0.3420 


0.0239 


0.5511 


1.5388 


5.3620 


2.7879 


0.6105 




v > v c (u) 


0.3420 


0.0239 





2.8591 


5.3620 


4.6651 


0.9171 




v < v c (u) 


0.2899 


0.0122 


0.5063 


0.0580 


8.6911 


< 1 




4 


v = v c (u) 


0.2899 


0.0122 


0.5063 


1.2051 


8.6911 


3.4427 


0.5717 




v > v c (u) 


0.2899 


0.0122 





2.0837 


8.6911 


8.4170 


0.9852 


3D #-chain (u = ug) 




v < v c (u e ) 


1/3 


1/3 


0.6180 


0.0613 


100/27 


< 1 




2 


v = v c (u g ) 


1/3 


1/3 


0.6180 


0.6180 


100/27 


1.8526 


0.4709 




v > v c (u s ) 


1/3 


1/3 





0.8229 


100/27 


2.4514 


0.6848 




v < v c (u g ) 


0.2071 


0.4307 


0.5511 


0.0211 


8.7231 


< 1 




3 


V = V c (ug) 


0.2071 


0.4307 


0.5511 


0.5773 


8.7231 


2.4203 


0.4081 




v > v c (u e ) 


0.2071 


0.4307 





1.2781 


8.7231 


3.7800 


0.6139 




v < v c (u e ) 


0.1918 


0.3393 


0.5063 


0.0180 


15.424 


< 1 




4 


v = v c (u e ) 


0.1918 


0.3393 


0.5063 


0.5758 


15.424 


3.5367 


0.4617 




v > v c (u e ) 


0.1918 


0.3393 





0.7984 


15.424 


5.9331 


0.6508 


3D globule (u > ug) 




v < v c (u) 





22- 1 / 3 


0.6180 





4 


< 1 




2 


v = v c (u) 





22- 1 / 3 


0.6180 


0.7637 


4 


2.4163 


0.6364 




v > v c (u) 





22- 1 / 3 





0.7637 


4 


2.4163 


0.6364 




v < v c (u) 





oo 


0.5511 





9.772 


< 1 




3 


v = v c (u) 





oo 


0.5511 


oo 


9.772 


2 


0.3041 




v > v c (u) 





oo 





oo 


9.772 


2 


0.3041 




v < v c (u) 





22-V3 


0.5063 





16 


< 1 




4 


v = v c (u) 





22-V3 


0.5063 


0.7637 


16 


5.8387 


0.6364 




v > v c (u) 





22-V3 





0.7637 


16 


5.8387 


0.6364 
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between them, both chains also have large parts that are not interconnected. Even for 
v > v c (u) fixed value (A*,B*) does not change, but then C* becomes equal to zero, 
meaning that the whole P 2 chain is covered with the P3 chain (which still has a lot of 
monomers in the bulk, far from the boundary in which P 2 lies). The regions v > v c (u) 
and v < v c (u) of the phase plane u — v, as well as the critical line v c (u), are additionally 
partitioned by the vertical line u = ug, so that each part obtained in such a way is 
characterized by different fixed point (A*, B*,C*,D*), corresponding to different phase. 
Coordinates of all fixed points are given in table Q], whereas in figure [3] one can see 
obtained phase diagrams for 6 = 2, 3, and 4 SG fractals. 

Shape of the line v c (u), as well as values of the exponent 0, show that the interplay 
between the intra- and inter-chain interactions in the system under study is quite subtle. 
In the 6 = 3 case, v c (u) decreases monotonically with u, meaning that stronger monomer- 
monomer attraction within the P3 chain eases its attaching to the P2 chain. However, 
for 6 = 2 and 4 fractals this is correct only for values of u up to ug, where v c (u) has 
its minimum. For larger values of u, function v c (u) monotonically increases with u, i.e. 
for u > ug intra-chain prevails inter-chain interaction and hinders attaching. Different 
behavior of the system on fractals with 6 = 2,4 and b = 3 is due to the peculiar topology 
of these lattices. Namely, although for u > ug polymer is in globular phase, compactness 
of that structure is not always the same. Only on b = 2 SG the globule is completely 
compact, i.e. its fractal dimension d9 is equal to the fractal dimension dj D of the 
lattice [21]. In the 6 = 3 and 4 cases d9 < d 3D , but this quasi-compactness is much more 
pronounced in the 6 = 3 case [22,23], which brings about different behavior of the system 
on the 6 = 3 SG fractal. Concerning the exponent 0, which can be taken as a measure 
of interconnection between the two chains, one can notice that it has different values on 
the critical line v c (u) and in the region v > v c (u), and in addition depends on intra-chain 
interaction parameter u (see tableland figure [3]). For each of the three studied fractals, 
in the range v > v c (u) the inequality 0(w < ug) > <p{u = ug) > (f>(u > ug) is satisfied. 
Such inequality could have been expected, since it means that when P 3 chain completely 
covers P2 chain, the number of contacts between them is smaller if structure of the P3 
chain is more compact. On the line v c (u), however, chain P 2 is only partially covered 
with P3, so that some similar conclusion is not plausible, which is indeed in accord with 
the calculated values of (p. Besides, it is interesting that for 6 = 2 and 4, the smallest 
value of is obtained for u = ug, which is not the case for 6 = 3 fractal. 

Finally, one should note that in the case 6 = 3, for the globular state of a solitary 
3D chain (u > ug), the coordinates of the corresponding fixed point are Aq = and 
Be = 00. Furthermore, a numerical analysis of function D^ r \ in the range v > v c (u) 
reveals that D* = 00. Nevertheless, the relation (M^) ~ \ r D and formula (12.111) are 
applicable, but with different meaning of A^. For the globule state of 6 = 3 fractal, it was 
demonstrated [22] that equations ( 1A.5I) and (1A.6jl in the vicinity of the corresponding 
fixed point (0, 00) have the following approximate form 

A (r+i) = 320(A (r) ) 3 (P (r) ) 6 , P (r+1) = 4308(A (r) ) 2 (P (r) ) 8 , (2.12) 
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Figure 3. Phase diagrams obtained for the model of two avoiding SAWs on 3D SG 
fractals with 6 = 2,3 and 4. The solid vertical line u = ug divides the u — v plane in two 
areas, corresponding to the phases in which the 3D SAW is either extended (u < ug) or 
collapsed (u > ug). Each of these two areas is additionally partitioned by the critical 
line v = v c (u). For v < v c (u) the two polymers are segregated one from another, and 
the system exists either as "3D extended SAW + 2D SAW" for u < ug, or "3D globule 
+ 2D SAW" for u > ug, whereas for u = ug precisely, 0-chain coexists with 2D SAW. 
For v > v c (u) the mean number (M) of contacts between the two SAWs scales with 
the mean length (N 3 ) of the 3D SAW as (A/) ~ (N 3 )'K Depending on the value of u, 
critical exponent <p has different values, which are presented within the corresponding 
areas. 
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Figure 4. The six restricted generating functions used in the description of all possible 
inter-chain configurations for the CSAWs model of the two-polymer system, within the 
r-th stage of 3D SG fractal structure. The 3D chain is depicted by green line, while 
the 2D surface-adhered chain is depicted by yellow line. 



from which it follows A, 



73+11 



9.772 and z/f 



In 3/ In 9.772 = 0.4819. Besides, 



vjyg — 2 

for v > v c (u), the inequality <C B^ r > is valid, so that RG equation (1A.8|) obtains 
the approximate form 



D (r+i) = 320A {r ' ) (5 (r) ) 6 C {r) ( J D (r) ) 2 . 

Then, from equations flUT} and (|21Sjl . follows (M^) = 2 D v (r) 9I £ ) 
that (M^) ~ X r D (for large r), with A D = 2. Finally, from (J2HJ, one obtains 
= In 2/ln 9.772 = 0.3041. 



(2.13) 

2(M^), implying 



3. The model of crossing walks 

In order to describe the physical situation when closer contact between the two polymers 
is possible, in this section we analyze the CSAWs model in which chains P2 and P3 can 
cross each other [14]. If we assume that chains interact only at the crossing sites, and, 
similarly as in the ASAWs case, introduce the weight factor w = e _ec//fcsT , where e c < 
is the energy of two monomers in contact, it turns out that the two chains cannot exist 
independently, even for extremely weak attraction (|e c | <C fc^T). Therefore, we define 
additional weight factor t = e~ et ' kBT , where tt > is the energy associated with two 
sites, visited by different SAWs, and both neighbouring a crossing site (see figure QJb)), 
so that unbinding transition can occur. To describe exactly all possible configurations 
of the two-chain polymer system, within this model we need to introduce nine restricted 
partition functions: A (r \ B (r \ C^, A { {\ A<£\ A ( £\ B[ r \ and B$ ] . Functions A^ T \ 
B^ and C^ r \ which correspond to one-polymer configurations are the same as in the 
ASAWs model (see figure [21 and RG relations (12.21) and (12.31) ). whereas the remaining 
six functions, which describe the inter-chain configurations, are depicted in figure HI and 
they are defined as 

4 W = E A W (iV2, N 3 , L, M, K)x^u L w M t K , i = 1, 2, 3, 4 , 

N 2 ,N 3 ,L,M,K 
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flW = £ Bf 3 (iV 2 , iV 3 , L, M, K)3$*x$*u L w M 1* , z = 1, 2 

N 2 ,N 3 ,L,M,K 



where .4] and i3 4 are the numbers of particular two-polymer configurations on the r- 
th fractal structure. For instance, A4 (N 2 , N 3} L, M, K) is the number of configurations 
in which the AVstep P 3 chain (with L intra-chain contacts) and iVVstep P 2 chain (with 
different entering end exiting vertices from P 3 chain) cross M times and have K pairs 
of sites belonging to different chains and neighboring the crossing sites. Functions A[ 
and satisfy the following recursion relations 

K = E ai W A N *B N *C N ° n A^ ft < Sfe , i = h 2, 3, 4 , (3.1) 

A/" j=l fc=l 

4 2 

S| = J] 6i(jV) A Na B Nb C Nc f[ A^ Aj n 5f Sfc , 2 = 1,2, (3.2) 

AT j=l k=l 

where M denotes the set of numbers Af = {Na, N b , N c , Na ± , Na 2 , Na 3 Na 4 , N Bl , N B2 }, 
and where we have used the prime symbol as a superscript for (r + l)-th restricted 
partition functions and no indices for the r-th order partition functions. The above set 
of relations (I3.ip - fl3.2l) . together with the previously introduced relations fl2.2p - fl2.4p for 
the functions A, B, and C, can be considered as the system of RG equations for the 
problem under study, with the initial conditions 

A^=x 3 , B^=xlu\ C^ = x 2 , 

A^ = x 3 x 2 w 2 , = A^ = x 3 x 2 wt , A^ = x 3 x 2 , (3.3) 

B[ 0) = B^ = xlx 2 w 2 u\ 

corresponding to the unit tetrahedron. Because the number of all possible configurations 
is extremely large, we have been able to find explicit form of the RG equations (13. ip - 



( 13. 2p only for 6 = 2 and 6 = 3 SG fractals (see Appendix B). For both cases numerical 
analysis shows that, for each considered value of t, there is a critical line w c (u,t) 
dividing the u — w plane into regions where the two polymers are either segregated 
(w < w c (u,t)) or entangled (w > w c (u,t)). Depending on the value of the intra- 
chain interaction parameter u, the area w < w c (u,t) is further partitioned into smaller 
regions, corresponding to various phases of the system (see figure [5]). To each of these 
area different fixed point of the general type 

(A*, B*,C*, Al, Al A3, At, Bl, B*) , (3.4) 

pertains. We describe general features of the fixed points and the corresponding phases 
in the three following subsections. 

3.1. Weak self- attraction of the 3D chain 



For each value of < t < 1, and small values of the interaction parameter 1 < u < uq, 
there is some critical value w = w c (u, t) such that 
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Figure 5. Phase diagrams in the space of interaction parameters for CSAWs model 
in the case of 6 = 2 and 6 = 3 SG fractal, for t = 0.5. In both cases the critical line 
w — w c (u,t) separates the u — w plane into the area w > w c (u,t) of entangled phase 
and area w < w c (u,t), in which the two chains are segregated. The latter area is 
divided by vertical line u = ug into regions, corresponding to three segregated phases: 
(i) 2D chain (always extended) and extended 3D chain (u < ug), (ii) 2D chain and 3D 
(9-chain (u = ug), and (iii) 2D chain and 3D globule (u > ug). One should observe that 
there appears the multi-critical point (full red circle) at the crossing of the 0-line and 
the critical line w = W c (u, t). For other values of t (0 < t < 1), the critical line w c (u, t) 
also monotonically decreases, for both 6=2 and 6 = 3 fractals. 



• For w < w c (u, t) the fixed point of the form 

{A E ,B E ,C*, 0,0, 0,^,0,0), (3.5) 

is reached. This point corresponds to the phase in which 2D chain and extended 
3D chain are segregated, since as it is approached, after some number r ^> 1 of 
RG iterations, the average number of contacts between the two chains, quickly 
becomes constant. Values of A E and B E are fixed values of the RG parameters for 
the solitary extended chain on 3D SG fractal, and they are presented in table [2} 
together with the values of C*, corresponding to 2D chain, which can exist only in 
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extended state. RG fixed point value A\ is equal to 0.1165 and 0.0779, for b = 2 
and 3 respectively, and they coincide with the values of D* for v < v c {u < u e ) case 
in the ASAWs model (see tabled]). 

For w = w c (u,t) one obtains the symmetrical fixed point 

(A E , B e , C*, A e C\ A e C\ A e C\ A e C\ B e C\ B e C*) , (3.6) 

which appears to be a tricritical fixed point. As one approaches this fixed point, the 
average number of contacts (M^) becomes infinitely large (although large parts 
of P2 and P3 are not in contact), and it turns out that it scales with the average 
length (N 3 ^) of the 3D chain, according to the power law 

(AfW) ~ (iV 3 (r) ) v . (3.7) 

To calculate the contact critical exponent (p, within the CSAWs model, we find 
the average number of contacts between two chains at the rth stage of fractal 
construction, through the formula 

V M /V 4 jWj-V 2 R( r A ~.N 2 N 3 L QI ,M f K 



(M (r) ) 



dAt ] "dBl Y 



1=1 

where Xi = A4 (i = 1, 2, 3, 4), X 5 = B±, and A 6 = B 2 . Since 

dx (r+1) « <9xj r+1) dxP 



-E *m ^ » i = l,...,6, (3.9) 
^ £1 <9X r) <9u> 



one expects, for large r, that d ^ behaves as AL where A v is the largest relevant 
solution of the eigenvalue equation 

«9x (r+1) ' 



det 



5x (r) J K S ij 



0, (3.10) 



where the asterisk means that the derivatives should be taken at the tricritical fixed 
point. From here follows (M^) ~ \ r , which together with (ivj ) ~ \ r U3 (where 
X U3 , as before, is the largest eigenvalue of the linearized RG equations for the bulk 
parameters A and B), and (13.71) . gives 

p = s^- (3 - n) 

For larger values of the interaction parameter w > w c (u, t), the RG parameters flow 
towards the fixed point 

(0, 0, 0, A\ = C*, 0, 0, 0, 0, 0) , (3.12) 

which describes the entangled phase of the two polymers, in which P3 chain is 
completely attracted to P 2 chain. 
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Table 2. The CSAWs model fixed points corresponding to the critical values 
w = w c (u,t) of the attraction parameter between the 2D and 3D chains, when 
< t < 1, for all possible states of the 3D polymer, together with the values of 
the critical exponent if. 



b 


A* 


B* 


C* 


A{ 


A* A* 


At 


Bl 


B* 2 


f 


extended 3D chain (u < ug) 


2 


0.4294 


0.0499 


0.6180 


0.2654 


0.2654 0.2654 


0.2654 


0.0308 


0.0308 


0.5428 


3 


0.3420 


0.0239 


0.5511 


0.1884 


0.1884 0.1884 


0.1884 


0.0131 


0.0131 


0.4973 










3D 


0-chain (u = ug) 










2 


1/3 


1/3 


0.6180 


0.0510 





0.0613 


0.2365 


0.2362 


0.6714 


3 


0.2071 


0.4307 


0.5511 


0.0810 


0.0310 0.0250 


0.0270 


0.3130 


0.3150 


0.6226 










3D 


globule (u > ug) 










2 





0.3569 


0.6180 











0.2206 


0.2206 


0.6261 


3 





00 


0.5511 











oo 


00 


0.6073 



3.2. Critical self- attraction of the 3D chain 

For u = ug the solitary 3D chain is in the state of the #-chain, for which (A*,B*) = 
(Ag, Bg), whereas the two-polymer system can be in the following phases: 

• For w < w c (ug,t) the corresponding fixed point is of the form 

(Ag, B e ,C*, 0,0, 0, Al,0,0) . (3.13) 

This is the case when the 3D ^-chain is segregated from the 2D chain chain. Fixed 
point values of A\ are: 0.0613 for 6 = 2, and 0.0211 for 6 = 3 fractal, equal to D* 
for the corresponding cases v < v c (ug) of the ASAWs model. 

• When w = w c (ug,t), the RG parameters tend to the fixed point 

(A e , Bg, C*,A\, A* 2 , Al, Al, B*, B* 2 ) , (3.14) 

which corresponds to the phase in which chains are not segregated anymore, but 
they are not yet completely entangled. In contrast to the w = w c {u < ug,t) case, 
for which symmetrical fixed point is obtained, values of Ai (i = 1, 2, 3, 4), as well as 
Bi and B 2 , are not mutually equal (Ai ^ AgC* , Bi ^ BgC*). The scaling relation 
(MM) ~ (N^ r) Y is satisfied, with <p given by (I3TTT1) . 

• For w > w c (ug, t) the RG parameters flow towards the entangled fixed point (13. 12ft . 

3.3. Strong self- attraction of the 3D chain 

When self-attraction of the 3D polymer is strong (u > u$), depending on the values of 
inter-chain interaction parameters, the following phases are possible: 

• For w < w c (u,t) the chains are segregated. Due to the large compactness of the 
3D chain, with (A*, B*) = (0, B G ), none of the configurations A\, A 2 , A 3 , A 4 , B\, B 2 
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can be accomplished, and the corresponding fixed point is 

(0,^,^,0,0,0,0,0,0). (3.15) 
The chains are completely separated. 

When attraction between the chains is critical, w = w c (u, t), the chains are partially 
entangled, and the fixed point 

(0, B G , G\ 0, 0, 0, 0, B G C\ B G C*) , (3.16) 

is attained. In this case the interaction between chains is sufficiently strong to 
connect them, but not strong enough to destroy the compactness of the 3D globule, 
so that all A* = 0. Again, the scaling relation (M^) ~ (N^)^ is satisfied for 
both 6 = 2 and 3, with (p given by (13.111) . However, while in the case 6 = 2 the 
coordinates of corresponding fixed point have definite values (and \ v can be directly 
calculated from linearizied RG equations for and Bj), in the 6 = 3 case some 
fixed point coordinates diverge, and calculation of A^ requires an additional effort. 
To be more specific, a numerical analysis of RG equations (13.11) and (13. 2ft reveals 
that A\ r) « A^C* -> 0, Bf ] « B^C* -> oo. In this situation the appropriate 
eigenvalue A^ can also be calculated, using the following transformation. If we write 
the relation (I3.9P in the form 

1 Efl^)^, (-) 



X\ r+l) dw fr{ \xt +1) dxf ] J Xf ] dw 
it can be shown that, when we keep only dominant terms in the RG equations, and 

dX (r+1) ( X {r) dX (r+1) \* 

in the derivatives — S^j-j then, the matrix elements — [7+ry — l tft °f the new 

OX \ X oX 1 

J \ 1 3 / 

eigenvalue problem are either equal to zero or to some finite constants (depending 
on C*), from which we find A^ = 3.9919, and therefrom (p = lnA^/lnA^ = 0.6073. 

• Strong inter-chain attraction w > w c (u,t) destroys the globule and completely 
attaches the 3D chain to the 2D chain. This entangled phase is again characterized 
by the fixed point (13.121) . 

In table [2] we have presented the numerical results for the crossover fixed points 
and the corresponding values of the contact exponent ip, obtained for the unbinding 
transitions from entangled two-polymer phase to segregated phases of 2D and 3D chains 
on the 6 = 2 and 6 = 3 3D SG fractals. These values are correct for all studied cases 
of t in the interval (0, 1). Varying the parameter t in this interval changes only the 
particular values of w c (u, t), but, for both 6 = 2 and 6 = 3, the function w c (u, t) for 
fixed t is monotonically decreasing function (see figure EJ). Dependence of w c (u,t) on 
t, when u is fixed, is presented in figure El for several values of u. As one can see, 
the limiting values t — and t = 1 are also included in this figure. However, in these 
cases different fixed points, from those obtained for < t < 1, can be reached, which is 
expounded in the following paragraphs. 

First, we analyze the value t = 0, which represents the limiting case, within the 
CSAWs model, when the energy e t (corresponding to the repelling of two different chain 
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Figure 6. Critical value of the inter-chain interaction parameter w c (u,t), depicted as 
a function of t, for three different values of intra-chain interaction parameter u, in the 
cases of b = 2 and b = 3 3D SG fractals. 

monomers, placed at sites which are nearest neighbours to a crossing site) is infinitely 
large. Starting with the initial values (13. 3p . it can be shown that, in the case of the 
6 = 3 fractal, the same fixed points of the RG equations (13. ip and (13. 2p . as for < t < 1 
are reached. However, for the 6 = 2 fractal, it can be seen, from the explicit form 
of the RG equations (jBJ^ - (jBT6l) . that t = leads to = = 0, for every r, 
x 2 , X3, u and w. This is due to the topology of this fractal, and a consequence is 
that the fixed point (Ae, Be, C* , C* , 0, 0, A\, A\, 0) corresponds to the critical values 
w = w c (l < u < uq, t = 0). The coordinates of this fixed point A E = A*, B E = B* and 
C* are given in the part "extended 3D chain (u < ug)" of the table EJ while Al = 0.1164, 
and the concomitant critical exponent is ip = 0.8439. The remaining fixed points are 
the same as for < t < 1. 

The second limiting value (t = 1) corresponds to the case e t = (when there is 
no repelling interaction). In this case, for all u, the critical value of the interaction 
parameter w is equal to w c (u,t = 1) = 1. This means that the chains can not 
be segregated, even for extremely small attraction between them. For both fractals 
6 = 2 and 6 = 3, the fixed points that pertain to the critical value w c , for u 7^ ug, 
are the same as for < t < 1. For u = ug the symmetrical fixed point is reached, 
(A e ,Bg,C*,AgC*,AgC*,A e C*,A e C*,BgC*,B e C*), in contrast to the case < t < 1. 
Values of Ag = A*, Bg = B* and C* can be found in the middle part of the 
table |2l while the values of the contact critical exponents are ip(b = 2) = 0.6102, 
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and <p(b = 3) = 0.5907. 

Finally, one should mention that recently, using the Monte Carlo renormalization 
group (MCRG) method, the contact exponent <p was calculated for b up to 40, for 
the case when the intra-chain interactions within the 3D chain are negligible, u — > 
[14]. Comparing the reported MCRG data ip MC (b = 2) = 0.5440 ± 0.0056 and 
v MC (b = 3) = 0.4969 ± 0.0024 with our exact findings <p(b = 2) = 0.5428 and 
ip(b = 3) = 0.4973, we can see that MCRG data are in excellent agreement with 
our exact findings. In [14] it was also demonstrated that ip MC , as a function of the 
scaling parameter b, continues decreasing with increasing b, and, it seems that in the 
fractal-to-Euclidean crossover region (i.e. in the limit b — > oo) it goes to the proposed 
zero Euclidean value. The inequality tp(b = 2) > (p(b = 3) is satisfied not only for 
weak intra-chain interactions [u < ug), but also for u > ug, as can be seen in table [2j 
Unfortunately, in the range u > ug, the MCRG calculation of ip is not feasible, so that 
prediction for the large b behavior of <p(u > ug), only on the bases of our exact data, is 
not possible. 

4. Summary and conclusion 

In this paper we have studied a system of two interacting chemically different polymer 
chains in a poor solvent. Such a situation can be modelled by two avoiding self-avoiding 
walks (ASAWs), as well as by two crossing self-avoiding walks (CSAWs). We assume 
that polymers are situated in fractal containers modelled by members of 3D SG fractal 
family, which are labelled by an integer b (2 < b < oo). We adopt that the first polymer 
(P3) is a floating chain in the bulk of 3D SG fractal, while the second (P2) is stuck to 
one of the four boundaries of the 3D SG fractal, which appears to be a 2D SG fractal. 
To take into account the intra-chain interaction of P3 polymer we have introduced 
the interaction parameter u = e _eu//fcsT , where e u < is the energy corresponding to 
interaction between two nonconsecutive neighboring monomers within the chain. In the 
case of ASAWs model, the two SAW paths cannot intersect each other, and we assume 
that two polymers interact when they approach a distance equal to a lattice constant. 
We associate the weight factor v = e~ ev ^ kBT with each such contact, where e v < is 
the appropriate energy of inter-chain interaction. On the other hand, in the case of 
CSAWs model, in order to describe the inter-chain interactions of P3 and P2, we have 
introduced the parameters w = e~ £c//fcsT and t = e~ et//feflT , where e c < is the energy 
corresponding to each crossing of SAWs, while e t > is the energy associated with a 
pair of sites, visited by different SAWs, which are nearest neighbors to a crossing site. 

To obtain the phase diagrams and the contact critical exponents between the two 
polymers, we have applied an exact RG method for the 3D SG fractals labelled by 
6 = 2,3 and 4, in the case of ASAWs model, and for fractals b = 2 and 3, in the case of 
CSAWs model. In both models, for various values of intra-chain interaction parameter 
u, a solitary 3D floating polymer chain can be found in one of the three possible phases 
(extended, 6*-phase, or globule phase), whereas a solitary 2D chain is always extended. 
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Depending on the values of the inter-chain interaction parameters (v in the case of 
ASAWs model, and w and t in the case of CSAWs model), the system can be either 
in the segregated phase, when the chains can be considered as almost independent, or 
in phases in which the number of contacts between the chains is comparable with their 
length (entangled phases). For both models, there is a critical line in the plane of the 
interaction parameters (v c (u) for ASAWs, and, w c (u, t), with fixed t, for CSAWs model), 
which divides it into areas corresponding to segregated and entangled phases. In the 
case of the ASAWs model, for v > v c (u), the average number (M) of contacts between 
the two polymers scales with the average length (N 3 ) of the 3D chain as (M) ~ (A3)*. 
Different values of the contact exponent <fi correspond to v — v c (u) and v > v c (u), and, 
in addition, <fi also depends on the strength of the intra-chain interaction parameter u. 
However, in all entangled phases large parts of the 3D chain remain in the bulk, beyond 
the scope of the inter-chain interaction, since the prohibition of crossings between two 
chains hinders their complete entanglement, even for extremely large values of v. On 
the contrary, in CSAWs model for w > w c (u, t) the two chains are completely entangled, 
while they only partially cover each other at the critical line w = w c (u,t), where the 
scaling relation (M) ~ (-^3)^ is satisfied, and where (p takes different values in the 
intra-chain interaction regions u < ug, u = ug, and u > ug. 

In the end, we would like to point out that for ASAWs model, in the space of 
interaction parameters, the arrangement of possible phases is approximately the same, 
as in the case of the surface-interacting polymer chain in a poor solvent in Euclidean 
spaces [15-17]. On the other hand, the obtained phase diagrams for CSAWs model, 
resemble the phase diagrams of the same surface-interacting chain problem, in fractal 
containers [23]. This similarity is not surprising, since in both models studied, one of the 
two interacting polymers is adhered to one of four fractal surfaces, and its monomers 
appear as a part of interacting surface (in the surface-interacting polymer problem). 
Here we may conclude that, our findings should be useful in making the corresponding 
3D models of the system of several interacting polymer chains in porous media. Besides, 
our results may serve inspiring in advancing theories of mutually interacting polymers, 
as well as for polymers interacting with boundary surfaces of homogeneous 3D lattices, 
in which case so far (to the best of our knowledge) an exact approach has not been yet 
made. 
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Appendix A. Renormalization group equations for the ASAWs model 

In this Appendix we give explicit RG equations for the model in which two chains avoid 
each other, for the cases b = 2, and b = 3 of 3D SG fractals. Equations for "bulk" 
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parameters A and B, as well as for the "surface" parameter C, were found in earlier 
works, and we give them here only for the sake of completeness. 
First, we give the RG equations for b = 2 3D SG fractal 



We note that first three equations were established for the first time in [19]. 
Next, we present RG equations for the 6 = 3 case 



A' = A 3 + QA 4 + 1QA 5 + 3AA 6 + 7QA 7 + 112A 8 + 112A 9 + QAA 10 + 8A 4 B + 3QA 5 B 
+ U0A 6 B + 292 A 7 B + A2AA 8 B + 332 A 9 B + 12 A 3 B 2 + 12 A 4 B 2 + 118 A 5 B 2 
+ 380A 6 B 2 + 80QA 7 B 2 + 6QAA 8 B 2 + 72A 4 B 3 + 352A 5 B 3 + 70AA 6 B 3 + 1728A 7 B 3 
+ 3AAA 4 B 4 + 15Q8A 5 B 4 + 8A8A 6 B 4 + 2QAA 4 B 5 + 3192A 5 B 5 + 320A 3 B 6 , (A.5) 

B' = A 6 + 12A 7 + AOA 8 + 60A 9 + 32A 10 + 28A 6 B + 88A 7 B + 22AA 8 B + 1Q0A 9 B 
+ A0A 6 B 2 + A9QA 7 B 2 + 59QA 8 B 2 + 17QA 5 B 3 + 7Q8A 6 B 3 + 105QA 7 B 3 + 88A 3 B 4 
+ 2QAA 5 B 4 + 253AA 6 B 4 + 1152A 4 B 5 + 1888A 5 B 5 

+ 5808A 4 B e + 1936A 3 B 7 + A308A 2 B 8 , (A.6) 
C' = C 3 + 3C 4 + C 5 + 2C 6 , (A.7) 
D' = 2A 6 D 3 + AA 7 D 3 + AA 6 D 4 + 2A 5 D 5 + 28A 6 D 3 B + AA 4 D 4 B + UA 5 D 4 B 
+ AA 4 D 5 B + 8A 4 D 3 B 2 + 5QA 5 D 3 B 2 + AAA 4 D 4 B 2 + AA 3 D 5 B 2 
+ 1AAA 4 D 3 B 3 + 3QA 3 D 4 B 3 + 12A 2 D 5 B 3 + 72A 3 D 3 B 4 + 132A 2 D 4 B 4 
+ 26AA 2 D 3 B 5 + 12A 6 D 2 C + 18A 7 D 2 C + 2A 4 D 3 C + 8A 5 D 3 C + 8A 6 D 3 C 
+ AA 5 D 4 C + 2A 4 D 5 C + 16A 4 D 2 BC + A8A 5 D 2 BC + 6AA 6 D 2 BC 
+ AA 2 D 3 BC + 8A 4 D 3 BC + A8A 5 D 3 BC + 4A 3 D 4 J BC + 8A 4 D 4 BC 
+ 8A 3 D 5 BC + \2A 2 D 2 B 2 C + 36A 4 L> 2 .B 2 C + 1624 5 Z> 2 B 2 C + 2AA 3 D 3 B 2 C 
+ 28A 4 D 3 B 2 C + 44A 3 D 4 fi 2 C + 8A 2 D 5 B 2 C + 96A 3 D 2 .B 3 C + QAA 4 D 2 B 3 C 
+ 152A 3 £> 3 .B 3 C + 24A 2 L> 4 J B 3 C + 2AAD 5 B 3 C + 512A 3 D 2 .B 4 C 
+ 88AD 4 B 4 C + 264A 2 D 2 J B 5 C + 320AD 2 J B 6 C + QA 4 DC 2 + 16A 5 DC 2 
+ 28A 6 DC 2 + 12A 7 DC 2 + 124 5 £> 2 C 2 + 18A 6 D 2 C 2 + 2A 3 D 3 C 2 
+ 6A 4 D 3 C 2 + 8A 5 D 3 C 2 + AA 4 DBC 2 + 6AA 5 DBC 2 + 56A 6 DBC 2 
+ \2A 3 D 2 BC 2 + SOA^EC 2 + SGA^SC 2 + 2AD 3 BC 2 + 6A 3 D 3 BC 2 
+ 2AA 4 D 3 BC 2 + lOOA 4 /^ 2 ^ 2 + SSA 5 /^ 2 ^ 2 + 6AD 2 B 2 C 2 
+ 18A 3 D 2 B 2 C 2 + 56A 4 D 2 B 2 C 2 + 4A 2 J D 3 5 2 C 2 + 20A 3 D 3 5 2 C 2 
+ 1Q0A 4 DB 3 C 2 + 3QA 2 D 2 B 3 C 2 + 32A 2 D 3 B 3 C 2 + 25QA 3 DB 4 C 2 
+ \32A 2 D 2 B 4 C 2 + 44AD 3 J B 4 C 2 + A 3 C 3 + 6A 4 C 3 + 10A 5 C 3 + 10A 6 C 3 
+ 6A 7 C 3 + 8A 3 DC 3 + 10A 4 DC 3 + 1AA 5 DC 3 + 10A 6 DC 3 



a 

D' 



A' 



B' 



A 2 + 2 A 3 + 2 A 4 + AA 3 B + 6 A 2 B 2 , 
A 4 + A A 3 B + 22 B 4 , 
C 2 + C 3 , 

2 D 3 5 + 6 £> 2 5 2 + 2 A 2 D C + A 2 C 2 + A D C 2 . 



(A.l) 
(A.2) 
(A.3) 
(A.4) 
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+ AA*D 2 C 3 + 6A 5 D 2 C 3 + 2A 3 D 3 C 3 + QA^D 3 C 3 + 8A A BC 3 + 16A 5 BC 3 
+ 20A 6 5C 3 + 28A A DBC 3 + 3QA 5 DBC 3 + 8A 3 D 2 BC 3 + 12A 4 D 2 BC 3 
+ AA 2 D 3 BC 3 + 1QA 3 D 3 BC 3 + 12A 3 B 2 C 3 + 18A 5 B 2 C 3 + 12A 2 DB 2 C 3 
+ AAA A DB 2 C 3 + Y2A 3 D 2 B 2 C 3 + 20A 2 D 3 B 2 C 3 + 2AA 4 B 3 C 3 
+ A8A 3 DB 3 C 3 + 2AA 2 D 2 B 3 C 3 + 2AAD 3 B 3 C 3 + 3A 3 C 4 + 10A 4 C 4 
+ 10A 5 C 4 + 4A 6 C 4 + A 2 DC 4 + 2A 3 J DC 4 + 2A 4 DC 4 + 12A*BC 4 
+ 8A 5 BC* + 4A 3 D J BC 4 + 18A 3 5 2 C 4 + 6A 2 DB 2 C* + 2A 2 DC 5 
+ AA 3 DC 5 + AA 4 DC 5 + 8A 3 DBC 5 + 12A 2 DB 2 C 5 . (A.8) 
Equations flA.51) and flA.6h were found in [22], and (IA.7I) in [20]. 



For the 6 = 4 case, equations are too cumbersome to be quoted here, and, they are 
available upon request to the authors. 

Appendix B. Renormalization group equations for the CSAWs model 

It can be shown, via direct computer enumeration of the corresponding paths within 
the generator of the 6 = 2 3D SG fractal, that RG parameters A±, A 2 , A3, A&, B\, and 
B 2 fulfil the following recursion relations 

A\ = A\ + A\ + AA\ + 2A 2 A 2 A 3 + AA\ + 2AA X A\ + 2^742^3^! + 2AA 2 A 3 B 2 

+ AA 2 A 1 B 2 + AA 2 B l B 2 + AAA 1 B 1 B 2 + 2A 2 B 2 + 2AA 1 B 2 2 + A\C + AA\C , (B.l) 

A' 2 = AA X A 2 + A\ + A 2 A 1 A 3 + AA 2 A\ + A 2 A 3 A± + AA 2 A\ + 2AA 2 A 4j B + AA 3 A A C 
+ A 2 A 3j Bx + AA 1 A 3 B 1 + 2AA 2 BB 1 + AA 2 B\ + A 2 A 3j B 2 + 4^^552 + A 2 A 3 C 
+ 2AA 2 5 1 E 2 + 3AA 2 B 2 + AA 2 C + AiA 2 C + AA X A 3 B 2 + 2AA 2 A 4 E 2 , (B.2) 

A 3 = A 2 A 1 A 2 + AA 1 A 3 + AA\ A 3 + AA 2 2 A 3 + A 2 A 2 A A + 2A 3 5 + AAA 3 BB 2 + AA 3 C 2 
+ 2AA 3 A±B + A 2 A 2 B 1 + AAiA 2j Bi + 2AA 3 BB X + 2^3^551 + AA 2 A 4 C 
+ A 2 A 2 B 2 + AA 3 A A BB 2 + A 2 A 2 C + A4 3 C + AA X A 3 C + AA X A 2 B 2 , (B.3) 

A' 4 = 2A 2 A 2 A 3 + 2 A4 2 A 4 + 2AA\B + 2AA 2 S + 2A\B + 6A 2 5 2 + 2A 2 3 BB 1 

+ 2AA\B 2 + 4A 2 S5 2 + 2AA 2 A 3 C + 2A 2 A 4 C + A 2 C 2 + AA 4 C 2 , (B.4) 

B[ = A 2 A\ + AA\A A + A4 2 5 + AA 2 S + A 2 3 A A B + 2,4 2 ,4 lJ B 1 + AA\B X + 8BB| 

+ G^E 2 + 2BB 3 + 2AA 2 5 2 + 8B 2 B l B 2 + ABB\B 2 + S^E 2 + 8BB X B 2 2 , (B.5) 

5 2 = A 2 A 2 A 3 + AA 1 A 2 A 3 + AA\B + AA 2 £ + A 2 A 4 5 + AA\B X + 10BB 1 B 2 + Gfifi 3 

+ 2A 2 A 1 5 2 + AA 2 5 2 + AA\B 2 + 12B 2 B 1 B 2 + 6BB 2 B 2 + lOfi 2 ^ 2 . (B.6) 

One can check, by inserting Ai = A 2 = A 3 = B\ = B 2 = and A4 = D, into equation 
( IB.4p for the function A 4 , that RG equation (1A.4P for the function D in the case of the 
ASAWs model is recovered. This is not surprising, since it follows from the definitions 
of A4 and D, and it is certainly also correct for the 6 = 3 fractal equations. However, 
here we do not quote the 6 = 3 RG equations because they are extremely intricate. 
For instance, equation for the parameter A\ has 2753 terms, and it is similar for the 
remaining Ai and Bi equations. 
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